Anyons and transmutation of statistics via vacuum induced Berry phase 
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We show that bosonic fields may present anyonic behavior when interacting with a fermion in 
a Jaynes-Cummings-hke model. The proposal is accomplished via the interaction of a two-level 
system with two quantized modes of a harmonic oscillator; under suitable conditions, the system 
acquires a fractional geometric phase. A crucial role is played by the entanglement of the system 
eigenstates, which provides a two-dimensional confinement in the effective evolution of the system, 
leading to the anyonic behavior. For a particular choice of parameters, we show that it is possible 
to transmute the statistics of the system continually from fermions to bosons. We also present an 
experimental proposal, in an ion-trap setup, in which fractional statistical features can be generated, 
controlled, and measured. 
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Anyons are quasiparticles that exhibit fractional quan- 
tum statistics which arises when bosons or fermions are 
confined in a two-dimensional space P, Q . Such excita- 
tions play a fundamental role in the physics of the frac- 
tional quantum Hall effect |3, 0] ■ Beside the theoretical 
applications, anyonic systems are promising for the im- 
plementation of fault-tolerant quantum computing 5]. 

In this paper, we show how typical bosons, specifically 
eigenstates of harmonic oscillators (number states), be- 
have as fractional spin particles, so-called anyons, when 
coupled to a two-level system (a fermion) through a 
Jaynes-Cummings-likeinteraction. 

We can interpret the origin of such an exotic behav- 
ior as a consequence of the entanglement capability of 
the interaction Hamiltonian (the entangled nature of the 
Hamiltonian eigenstates). When the coupled system is 
subjected to an adiabatic evolution that mixes the or- 
thogonal modes of the harmonic oscillator, the entangle- 
ment provided by the Hamiltonian can be regarded as 
a two-dimensional constraint on one of the two subsys- 
tems, allowing for bosonic excitations of fields (photons 
or phonons) to behave like anyons. This fractional be- 
havior is clearly manifested in the geometric phase |^ 
acquired by any eigenstate of the coupled system under 
a cyclic adiabatic evolution, and it is possible to relate 
the geometrical phase with the statistical factor Q . 

The bosonic field behaves, in an analogous way, as a 
spin-1 particle when there is no interaction (in the sense 
that the geometric phase acquired during the evolution 
is the same as that acquired by a spin-1 particle in a 
magnetic field), to a spin- 1/4 particle when the two sys- 
tems are maximally entangled, going through all the in- 
termediate spins depending on the specific degree of en- 
tanglement between the two subsystems. More specifi- 
cally, it is possible to simulate anyons with 
{m — 1, 2, ...) statistics or even to transmute continually 
the statistics of the system from Fermi to Bose. 

We also present a proposal, employing trapped ions, in 



which fractional statistical phases can be generated, ma- 
nipulated, and tested. This proposal, in addition to its 
experimental interest, provides a clearer physical frame- 
work to interpret the presented theoretical results. 

Initially, we consider an effective two-level (fermionic) 
system coupled nonlinearly to two quantized bosonic 
fields, through the so-called m-quantum Jaynes- 
Cummings model jH. In the rotating- wave approxima- 
tion, this Hamiltonian is given by (h — 1) 

H — i/a^ a+iyb^b+—crz + X,n cr+ (a)™ + ct- (a''^)™ , (1) 

where a, = |t) (Th li) (i| , a+ = |T) (i| , and (7_ = |i) (T| 
are the usual Pauh pseudospin operators (|t) and ||) are 
the excited and ground states of the two-level system, 
respectively). Am is the effective m nonlinear coupling 
constant, (a) and 6^ (6) are the creation (annihilation) 
operators of the bosonic fields with frequency v. 

In the interaction picture and in a rotating frame 
(through the unitary transformation exp [iA,„crzi/2]), H 
can be rewritten as 



H 



(T+ (a)™ + (T- (a^y 



(2) 



with Am = Lu — mj/ being the effective detuning between 
the bosonic and fermionic systems. Field mode a is or- 
thogonal to mode b and, initially, the two-level system 
interacts only with the former one. 

The eigenstates of the joint system, as- 
sociated to the eigenvalues ±A (A — 



(A,„)V4+(A„0'(n-^m)!/n! 

by 



are described 



where 



-(CT|T)|n),±QU)|n + m)J|n'),, (3) 
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and \n)i {I = a,b) is the Fock state of the bosonic field 

Let us consider the case in which the Hamiltonian 
described by Eq. ^ is slowly transformed into a 
parametrized one, H(6',(/)) = [/(6l, 0)HC/t(6', (/>). This 
transformation is achieved by externally driving the in- 
teracting systems, and if it is done adiabatically, i.e., 
slowly compared to the typical time scales of the prob- 
lem, the eigenstates of the initial Hamiltonian follow the 
transformation, as defined in the adiabatic theorem. In 
the end of a closed cycle (i.e., 9 and (j) are cyclically 
varied and then brought back to their original value), 
each eigenstate goes back to the original one, except for 
a phase factor of geometrical nature e^'^"-"' [ajlll- 

We are interested in the parametric transformation of 
Hamiltonian Q obtained through the unitary operation 



U (6*, (f>) = cxp{—i(j)Jz) exp{~i9Jy), 



(4) 



where J, 



(a'I'a — b''b) /2 and Jy — i {ab"^ — a'l'fe) /2 are 
the Schwinger angular momentum operators. In this 
case, for each value of the parameters 9 and 0, this trans- 
formation results in the two-level system interacting with 
a linear combination of the two bosonic fields, a and b 

A suitably slow variation of the parameters 9 and re- 
sults in the adiabatic evolution of the eigenstates \E'^ ^, ^ 

of the initial Hamiltonian H in Eq. (O . When 9 and 
are eventually brought back to their initial values. 



n,n' 



acquires a Berry phase, which is given by 



n,n' 



71, n' 



^in-n'+('^]^-H^^^-±^^], (5) 



(A2 + Ai 



where ft is the solid angle subtended by the cyclic path 



in the Poincare's sphere. Both eigenstates 



*„,„' > and 



n' I acquire the same geometrical phase. 
We are interested in the case of zero initial excitation 
in the bosonic fields (n = n' = 0), and the most interest- 
ing scenario is achieved when = 0, i.e., the resonant 
interaction. In this regime, the harmonic oscillator and 
the two-level particle are exchanging m excitations and 
the geometric phase reduces to 70,0 — x^- l^otice that, 
in this expression, each order of the non-linear interac- 
tion in ^ is contributing with a factor of 1/4 to the 
total geometric phase. The 1/4 factor originates in the 
fermion-boson interaction as shown in [lo|. while the m 
factor is typical for a collective behavior. In this sense, 
the Hamiltonian can be regarded as a multibody in- 
teraction, between the two-level system and m-bosonic 
particles. 



The closed path described by the transformation Q 
can be regarded as two consecutive exchanges of bosonic 
excitations in the two possible modes, which here play 
the role of two different spatial configurations of one any- 
onic particle. For example, initially the two-level sys- 
tem is exchanging m excitations with mode a. Then, for 
9 = 7:, the exchange involves only mode h, and finally for 
9 = 27r, a complete cycle is done and the interaction is 
back to mode a. This is analogous to the case of one elec- 
tron orbiting around a magnetic fiux tube in the original 
Wilczek work 0]. The analogy with the classical rota- 
tion in real space is clearer in the trapped ions example 
presented later, where modes a and b can be chosen as 
nothing but spatial vibration modes of the trap (see Ref. 
|l2l | for the implementation of fractional dynamics in the 
context of Bose-Einstein condensation) . 

In analogy with the physical two-dimensional space, 
we exploit in our system the parametric space of the 
Hamiltonian. When the parameters of the Hamiltonian 
H(0, 0) are changed and eventually returned to their orig- 
inal configuration, the wave function remains the same, 
except for the phase factor e*"^'^, where a = 7o,o/27r is 
called the statistical factor. In a 3-dimensional space, the 
rotation group satisfies a peculiar non-Abelian algebra, 
which allows only for the Bose (a = 0) or Fermi statis- 
tics (a = 1). On the other hand, if we restrict ourselves 
to 2-dimensional rotations, the correspondi ng g roup can 
generate a broader class of (braid) statistics [13] . It turns 
out, in fact, that the impossibility of using rotations in a 
third dimension reduces the number of constraints on the 
symmetry of the wave function, allowing for fractional 
values of a. 

We can interpret the eigenstate of the system described 
by Eq. Q (an entangled state of the two-level system 
and the field) as a kind of two-dimensional confinement, 
in the sense that only the rotation contributes for 
the Berry connection in Eq. ij^l, introducing the frac- 
tional behavior. To deepen our understanding of the role 
played by the entanglement of the eigenstates of the sys- 
tem (O, we will analyze the statistical factor a in dif- 
ferent regimes, ^ Am, A^ ~ A„i and A„j = 0. In 
general, the statistical factor is given by 



m Q. (Am) JTi! 
T2^(A2-hA%^ 



(6) 



and presents a dependence on the order m of the nonlin- 
ear interaction, on solid angle fi, and also on the detuning 

Am- 

Let us consider the far off-resonance case. Am ^ Am, 
which means that the harmonic oscillators and the two- 
level system are not exchanging energy, i.e., the eigen- 
state of the Hamiltonian is completely separable. This 
is the most trivial case, where a = and the excita- 
tions of the harmonic oscillators behave as bosons, as we 
should expect for a noninteracting system On the 
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FIG. 1: (a) The fractional part of the geometrical phase 
|7o,o (x^) ^} '^^'^ C-*) linear entropy of the reduced 
two-level system in the state (|HJ as function of the effective 
detuning Am- 



other hand, in the resonant case (A^ = 0) the eigen- 
state of Hamiltonian Q is maximally entangled and 
the statistical factor a = tI^- Therefore, the frac- 
tional features of the system's geometric phase depend 
crucially on the entangled form of the system's eigen- 
state. In Fig 1(a), we show the ratio 70,0 (t'St) ^ 
a function of the detuning A™, and in Fig. 1(b) we 
show the linear entropy of the reduced two-level system 
in state l|2Jl Sf = 1 — Tr/ {Tr^pf' (where / and b 
stand for fcrmionic and bosonic variables, respectively, 
and p = l^ao) (^o^o|) ^ function of the detuning 
Am- Comparing these two figures, we can see that, when 
the detuning increases, the amount of entanglement de- 
creases and consequently the fractional features of the 
system also decrease, until they disappear completely for 
separable eigenstates. 

It is interesting to note that, if we vary 9, (j) in a, way 
that the cyclic path in Poincare's sphere encloses the 



whole sphere (fl = An), the statistical factor, for Am = 0, 
turns out to be a = ^ , if we choose a path that encloses 
2/3 of the sphere, a = ^; and so on. In this way, we can 
simulate anyons with statistics. Let us con- 

sider the interesting particular situation (for A^ ^ 0) in 
which we choose m = 2 and vary the parameters 9, (p 
in order to obtain Q, = 47r. In that way, the statistical 
factor depends only on the detuning A2 . In the resonant 
case (A2 = 0), we have a — 1 and the system obeys the 
Fermi statistics. On the other hand, in the far detuned 
case (A2 » 1) the statistical factor tends to zero {a 0) 
and the system obeys the Bose statistics. In this particu- 
lar situation, varying the detuning A2 we can transmute 
continually the statistics of the system from Fermi to 
Bose, going through the anyonic statistics. The statisti- 
cal factor a depends on A2 with the same pattern showed 
in Fig. 1(a) for m = 2 (dashed line). We note that for 
the case m — 1 the value a = I cannot be reached. In 
this latter case, the statistical factor can vary continually 
from to 1/2 depending on the detuning Ai for the fixed 
solid angle = 47r. 

The scenario presented here for the simulation of 
one anyon statistics can be generalized for two anyons. 
Since the system composed by a fermion interacting 
with two bosonic fields behaves like one anyon, the 
generalization to two anyons can be achieved by adding 
two more bosonic fields (c, d). Let us suppose that the 
dynamics of the system is governed by the parametrized 
Hamiltonian (in the interaction picture) H(0, 0) = 
XmU{9, 0) [a+ {ar (c)'" + a. (at)'" {c^r]UH9, 0), 
under the unitary transformation U {9, (p) — 



exp I 

jlk 
•J z 

tum 
and 
the 



and Ji^ are the 



(J- 



31% (where 
y Schwinger angular momen- 

operators for modes I and fc). Under cyclic 
suitably slow variation of parameters 9 and 0, 
eigenstates of this new Hamiltonian, given by 

(|T)|0)jO),±|i)HJm)J|0)jO),/V2, 



evolve adiabatically. If 9 and are eventually brought 
back to their initial values, ^0^^ acquires the Berry 

phase 7o = which corresponds to twice the phase 

obtained in the previous case. Consequently, for the 
two-anyons case we have the statistical factor a = ^ 57 ' 
which is exactly twice as much as the one described in 
Eq. ® for A™ = 0. 

Finally, we discuss how to implement, in a physical 
context, Hamiltonian and the unitary transformation 
(gl). To this end we consider one single ion in a two- 
dimensional harmonic electromagnetic trap in the x-y 
plane, with degenerate frequency v. The ion has two ef- 
fective electronic states, ||) and ||), separated by the fre- 
quency Wo and coupled by the interaction with an effec- 
tive laser plane wave propagating initially in the x direc- 
tion, with frequency and wavevector k l = {lul/c) Ic - 
In this configuration, only the ionic motion along the x 



4 



axis will be modified and the Hamiltonian of this system 
is given by [l^ 



H 



LOq 



^ik L - r -iiAjLt+iipL 



H.c. 



^ (7) 

where r is the vibration direction of the ion, g the effec- 
tive coupling constant for transition ||) < — > ||), and (pL 
is the phase of the laser The effective laser beam 

is tuned to the m-th red vibrational sideband of the ion, 
i.e., it is detuned hy S = luq — uj^ = mv from the |t) 
< — > 11) transition. In the interaction picture and in the 
frame rotating at the effective laser frequency yf,, the 
Hamiltonian of this system can be written as [ij, [ig 



This laser interacts with the ion for a time r ^ l/A™, 
during which its propagation direction is rotated and lpl 
is cyclically changed. The variation rate of these param- 
eters must be much smaller than the effective coupling 
Am to allow the adiabatic regime assumed in our ap- 
proach. When the laser is turned off, the system has 



evolved to state [e'''">-« |t) + U)) |0)^ |0)fc/\/2; here we 
have used the fact that both eigenstates 



) and 



H = /m(77,ata) e'^-a+ (a)" + H.c. 



(8) 



where /,„ (77, at a) = \ge '''/^ 



is the effective coupling, 77 = 



1^1=0 l\(l+m)\ I" ) \"') + 



H 



2Mv is 



the Lamb-Dicke parameter, and Ai is the ion mass. 
When we assume the so-called Lamb-Dicke regime 



77 ^ 1, we have /,„ (77, a^a) 



2(m!) 



9 (ill) 



/2 = A 



Different values of the nonlinear interaction order m can 
be reached by the choice of the m-th red vibrational 
sideband of the ion. In this way we can implement 
Hamiltonian If we consider the effective pumping 
laser beam tuned not exactly on the m-th vibrational 
sideband of the ion, we can obtain the off-resonance 
case of Eq. |^ with a small detuning limited by the 
frequencies of the neighboring side bands {Am << i^)- 
In this context, the statistical transmutation of the 
system can be investigated varying the solid angle 
that depends on the path choice for parameters 9, 4>. 
The unitary transformation can be implemented 
by (i) the physical rotation, in the x,y plane, of the 
propagation direction of the effective laser field, in a way 
that k L = {ujlIc) [cos (6'/ 2) -I- sin (6'/ 2) ~y] , together 
with (m) a suitable variation on phase ip^ (note that 
ifL = 2y)). Therefore, by the control of the propagation 
direction of the laser and its phase, it is possible to 
implement physically the unitary transformation Q). 

The fractional phase acquired by the ion due to uni- 
tary transformation can be measured by a Ramsey- 
type interferometer similar to the one suggested in 
Let us start with the assumption that the vibrational 
modes of the ion in both directions x^y are cooled to their 
ground states, and its electronic internal levels are pre- 
pared in a superposition (|t) -I- |0)^ |0)j/ This 
superposition can be generated by a carrier-type laser 
pulse [3, which corresponds to the choice 777 = (the 
laser is tuned in the jl) < — > ||) transition frequency, 
ldq = lul), and does not affect the vibrational modes of 
the ion. The next step consists of turning on a laser 
tuned to the m-th red vibrational sideband and initially 
aligned in the x direction with reference phase ipL = 0- 



acquire the same geometrical phase 70,0 Siiid 

T> |0>a |0>6 = (l*ao) + 1*0,0))/ V2- It is clear that the 
variation of parameters 6 and ipL must be designed such 
that the end of its cyclic path coincides with the end of 
the j-th Rabi cycle performed in the time r = 27rj/ Am 
U 3> 1). Finally, we turn on again a carrier type 
laser pulse (m = 0) with phase (p^ — 7r/2 during time 
c-i = tt/ {2g). After this rotation we perform a fluores- 
cence measurement of the electronic states of the ion 14] , 
with probability P| = (1 — cos (70,0)) /2 to measure level 
II). In this way, we can measure the fractional phase 
introduced by the above procedure. 

In this paper, we have presented a method to simulate 
the dynamics of anyons via the Jaynes-Cummings model. 
We have shown how to simulate anyons with 
statistics or even how to transmute continually the statis- 
tics of the system from Fermi to Bose, going through the 
anyonic statistics. Such fractional features depend cru- 
cially on the entangled form of the system's eigenstate 
which works as a two-dimensional confinement. We also 
provide a proposal for the physical implementation of the 
above ideas in the context of trapped ions phenomena. 
It introduces a novel possibility for the investigation of 
fractional statistics in a well-controlled way. 
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